Slowly rotating fluid balls of Petrov type D 



o 
o 

(N 
O 

Q 



> 
o 

(N 

o 

"o 

i < 
5h . 

bo 



Michael Bradley*, 1 Daniel Eriksson''', 1 Gyula Fodor-f, 2 and Istvan Racz§ 2 

'Department of Physics, limed University, SE-901 87 Umed, Sweden 
2 KFKI Research Institute for Particle and Nuclear Physics, H-1525, Budapest 114, P.O.B. 49, Hungary 

The second order perturbative field equations for slowly and rigidly rotating perfect fluid balls 
of Petrov type D are solved numerically. It is found that all the slowly and rigidly rotating per- 
fect fluid balls up to second order, irrespective of Petrov type, may be matched to a possibly 
non-asymptotically flat stationary axisymmetric vacuum exterior. The Petrov type D interior so- 
lutions are characterized by five integration constants, corresponding to density and pressure of 
the zeroth order configuration, the magnitude of the vorticity, one more second order constant and 
an independent spherically symmetric second order small perturbation of the central pressure. A 
four-dimensional subspace of this five-dimensional parameter space is identified for which the so- 
lutions can be matched to an asymptotically flat exterior vacuum region. Hence these solutions 
are completely determined by the spherical configuration and the magnitude of the vorticity. The 
physical properties like equations of state, shapes and speeds of sound are determined for a number 
of solutions. 
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I. INTRODUCTION 

In [l[ a second order formalism for slowly and rigidly 
rotatingstars was developed by Hartle. This formalism 
was in [2[ applied to rotating white dwarfs and neutron 
stars using the Harrison- Wheeler and Tsuruta Cameron 
V 7 equations of state and in Q to the case with con- 
stant energy density. In [i[ the second order formalism 
is compared with numerical solutions of the full Einstein 
equations. For a review of relativistic rotating stars see 
0. In 0] global models for slowly rotating bodies in the 
post-Minkowskian approximation are treated. In a re- 
cent paper 0] second order perturbation theory for the 
matching of general stationary axisymmetric bodies to 
an asymptotically flat vacuum has been put on a more 
solid mathematical ground and the exterior metric is de- 
termined to second order. 

In this paper we use the Hartle formalism to study 
perfect fluids of Petrov type D. This condition will be 
used instead of an equation of state. It was shown in |8[ 
that physically realistic rotating fluid balls cannot be of 
Petrov types II, III, N or so the only possible cases are 
of Petrov types D or I. Hence it is of interest to closer 
study the properties of Petrov type D solutions, being the 
only possible algebraically special solutions. Also, since 
in the non-rotating spherically symmetric case all inte- 
rior solutions are of Petrov type D or 0, one might hope 
to find physically interesting interior solutions of Petrov 
type D also in the axisymmetric case, at least for slow ro- 
tation. However, the quadrupole moment of the rotating 
configuration will typically deviate from that given by 
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the Kerr metric and hence its exterior metric cannot be 
Kerr Q . It is easily verified that such an exterior metric 
is not of Petrov type D. 

The field equations to second order in the small ro- 
tational parameter will be solved numerically using 
fourth order Runge-Kutta. The system reduces to a 
closed subsystem of six first order differential equations. 
There are also two more differential equations for two 
further dependent variables which do not appear in this 
closed subsystem. Assuming regularity at the centre the 
solutions of this closed subsystem depend on four con- 
stants of integration, corresponding to zeroth order cen- 
tral density and pressure, the magnitude of the angular 
velocity and one more second order small constant. Due 
to scaling invariances we need only consider a two dimen- 
sional subspace of the solution space. The solutions are 
then matched to a second order axisymmetric_vacuum 
solution using the Darmois-Israel procedure [Tfjj . This 
metric includes the general second order asymptotically 
flat stationary axisymmetric vacuum solution as a spe- 
cial case. The interior solutions that can be matched to 
this vacuum form a three dimensional subspace of the 
space of solutions. One more freely specifiable parame- 
ter, associated with an independent spherical symmetric 
second order small change of the central pressure, is ob- 
tained from the solution of the two remaining equations. 
Hence the rotating configuration for the asymptotically 
flat subclass is determined by the spherically symmetric 
configuration (including a possible second order change 
of the central pressure) and the magnitude of the angular 
velocity. 

The paper is organized as follows: In section [TT] the 
method is briefly described and the field equations are 
presented, along with the Petrov type D condition. Fi- 
nally, the second order vacuum metric is given. The 
matching procedure is described in section IIIII and the 
integration constants for the vacuum solution are solved 
for in terms of the values of the interior solution on the 
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matching surface. In section ITVl the equations are rewrit- 
ten in a form suitable for numerical integration. The re- 
sults of the numerical runs are given in section |V] First 
the program is checked against the exact Wahlquist so- 
lution and then the subset for which the solutions are 
asymptotically flat is determined. Properties like shape, 
equation of state and speed of sound are then determined 
for a number of solutions. 



possess the form 

« Q = (lA^rn, 0, 0, 0) = ((1 - h)/A, 0, 0, 0) 

which also implies that the shear of the fluid is zero so it 
rotates rigidly. 

A. The field equations 



II. PRELIMINARIES 

To second order the metric of a slowly rotating ax- 
isymmetric object, both in the interior fluid region and 
the outside vacuum region, can be written as 
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ds 2 = (l + 2h)A z dt z - (l + 2m)^rdr - 

(1 + 2k)r 2 \d6 2 + sin 2 6 (dip - Lodtf] , (1) 

where lo is first order and h, to are k are second order 
in the rotational parameter [l[ . The requirements of reg- 
ularity at the centre and asymptotic flatness imply that 
the first order function lo depends on r only. The second 
order functions h, to, k can be given as 



h = h + fi2P2{cos8) 
too + to 2 P2(cosi 



to 
k 



= fc 2 P2(cOSI 



(2) 



where /iOi m o and /12 , ^2 , are functions of r only and 
P2{x) = i(3a; 2 — 1) is the second order Legendre poly- 
nomial. This result follows from reflection symmetry in 
the equatorial plane, from that the equations for h, to 
and k separate with the ansatze h — Y^iLo hi(r) Pi (cos 9) 
etc., and from the fact that there are no inhomogeneous 
terms containing lo in the equations for hi , ki and rrii for 
i > 2. For more details see 

The matching of the two spacetime regions happens 
via the application of a coordinate transformation ip — > 
tp + Qt in the fluid region. In addition to this, we can also 
rescale the interior time coordinate first by a constant 
C4 while matching the spherical zeroth order solutions, 
and then later by a second order small constant when 
doing the matching of the corresponding rotating space- 
times. These yield then the coordinate transformation 
t — > Ci(\ + c^)t, The first of these coordinate transforma- 
tions says that the inner fluid region rotates with respect 
to the distant stationary observers with angular velocity 
f2. This parameter Q, is considered to be the small ex- 
pansion parameter with respect to which lo is first order 
and the other corrections h, to, k are second order. 

The matter content of the interior is modelled by a 
perfect fluid 

Tab = (p + p)u a u b - pg ab . 

The coordinate system used in (TTJ) is assumed to be co- 
moving with the fluid, i.e. the 4-velocity is assumed to 



The basic field equations relevant to various orders can 
be listed as follows . If no equation of state is specified 
then the only equation one gets to zeroth order of the 
rotational parameter is the pressure isotropy condition 
G 1 1 = G 2 2 which reads as 
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d 2 A d(rA) d(B/r) 



dr 2 



dr 



dr 



A 







(3) 



Making use of G° = T° and G 1 1 = T 1 1 the energy 
density and pressure of the non-rotating configuration 
reads as 

d(rB 2 )' 



1 



Po 



Pa 



1 



1 



dr 

B 2 d(rA 2 ) 
A 2 dr 



1 



(4) 
(5) 



To first order in the rotation parameter the only relation 
B du>\ 



follows from G 3 = 



d_ 

dr 



A dr J dr 







(6) 



The second order Einstein equations yield the following 
four conditions. From G 1 2 = one gets 



d 1 dA 
r— (h 2 + k 2 ) + r (h 2 - m 2 ) — — h 2 - m 2 = 



- -- Adr - ■- --P) 

The pressure isotropy condition in the angular directions, 
G 2 2 = G 3 3 , gives 

6 (h 2 + to 2 ) 



,B 2 



A 2 



duj 
dr 



+ 4r 3 w 2 



B d 



A dr 







(8) 

The equality of the pressure in the angular and radial di- 
rections, i.e. G 1 1 — G 2 2 gives two equations. After elim- 
inating the derivative of h% using ([7]) one obtains from the 
P2(cos6') part 



B 2 dA 
A dr 



dk 2 
dr 



m 2 



, d (B 
2r 2 Bh 2 — - 
dr \ r 



+ 



1 4 B 2 



3 A 2 

while the ^-independent part takes the form 
d 



Tr) =° 



(9) 



Qr A B 



dr 



12A 2 m 



l -A 2 B^ 

r dr 



3r*B 2 [ *L 



dr 



9 d(r 2 A 2 ) dm 
dr dr 

, , , d (B 
+ 4r 3 Lu 2 AB— — 



dr V A 



= . 
(10) 
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The energy density function can be decomposed as p = 
P0 + P2, where p 2 = P2O+P22-P2(cos0) and p 20 and p 2 2 are 
second order small functions of the coordinate r given as 



P20 



B 

6r 2 A 



ir 3 uj 2 



D 

dr \ A 



^4 d / \ 



B ( duj 



A V dr 



and 



P22 



2(3A 2 h 2 + r 2 tj 2 ) r 



1 - B 2 + r 2 — ( B — 



d 
dr 



.dB 
dr 



(11) 



The analogous decomposition of the pressure is defined 
by P = Po + P2 = Pa + P20 + P22P2(cos9), where 



zeroth order pressure and density, given by ([5]), one gets 
24rA 4 

dr z 

dAdB 



dr 2 

AA 2 (B 2 - 1) -4r 2 AB 



dh „ , . dA 
— - + 2ArAm — x 
dr dr 



4 2 A 2 B^ 
dr dr dr 2 



+ 



ir 3 A 5 



dA d (B 



A 2 



2r 2 



dr dr \ A 2 
,d 2 B 2 



dmo r 2 uj 2 d ( B 2 



dr 



B 2 dr V A 2 



dr 2 



9 t dA d ( B 2 \ , 9x 
2rA ^dr{-A 2 ) + ^- B ^ 



— J - 12 A 2 m Q 



= 



(19) 



B. The Petrov type of slowly rotating fluids 



B 2 



P20 



and 



6r 2 A 2 



P22 



. o dhn d , 

l2rA —— - 12m — (rA 



dr 



dr 
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2B 
ir~A 



(3^ 2 + rV)A( 



The existence of a barotropic equation of state p 
is equivalent to 



dp dp dp dp 
d6 dr d6 dr 







(13) 

(14) 
Pip) 

(15) 



which is a geometric condition ensuring the coincidence of 
the constant pressure and density surfaces. Substituting 
the decompositions p = po + P20 + P22-P2(cos #) and p = 
Po + P20 + P22^2( c °s 0) into this relation yields, up to 
second order, 



P22' 



dpo 
dr 



- P22 



dpo 
dr 



(16) 



which is identically satisfied in virtue of the above field 
equations. 

It seems to be plausible to require the equation of state 
to be independent of the angular velocity. This condition 
reads as 



P2 
P2 



dpo 
dp Q 



(17) 



The ^-dependent part of this relation is equivalent to 
(fl"6|) . while the spherically symmetric part gives the rela- 
tion 



P20- 



dpo 
dr 



P20- 



dpo 
dr 



(18) 



Then by the substitution of the expressions p 2 Q and p 2 o, 
given by (|11[) and (fTB")) , together with po and po for the 



The spherically symmetric field equation ^ is usually 
complemented by a choice of an equation of state for the 
fluid. Since spherically symmetric static spacetimes are 
always algebraically special, we do not assume any spe- 
cial equation of state for the non-rotating base solution. 
As we will see shortly, the deviation from algebraically 
special cases can arise first when considering the second 
order terms in the rotational parameter. Here we require 
the interior solution to remain Petrov type D for slow ro- 
tation, thereby completing the system of field equations 
©, ©, ©, ©, ©, (HHD and (HU) by a further condi- 
tion, which in some sense plays the role of an equation of 
state. 

In order to calculate the Petrov type we need a suitable 
null tetrad. Up to second order an orthonormal tetrad 
can be given as 



0A 1 _ 



I' 1 



2A 2 

(0,B(1 - m),0,0) 

0,0,-(l -k),0 
r 



ruj 1 
— sin 6», 0,0, — — 
A z r sin 



,0,0,0 



-1 



2A 2 



sin 



• + k 



From this we form the null tetrad by the relations 

\/2^ = eft + eg 
V2k" = eft -eg 
y/2 = e^ + ie^ . 

Then the components of the Weyl spinor are given as 

*o = C abcd k a m b k c m d * 3 = C abcd k a l b m c l d 
* 1 = C abcd k a l b k c m d ^4 = C abcd m a l b m c l d 
* 2 = C abcd k a m b fh c l d . 

Since ^1 — and ^ = (even in case of fast rotation) 
the Petrov type is determined by the multiplicities of the 
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roots of the algebraic equation for the complex number 



*o + 6* 2 a 2 + *4a 4 = . 



(20) 



The Petrov type is D if there are two double roots, i.e. 



9*, = ^n* 



0*4 • 



(21) 



We note that the Petrov type can also be D if ^ = ^4 = 
and $2^0 but since then the equation of state can be 
shown to be p — —p (see Q) we only deal here with the 
more general case (|2"Tj) . Considering the other possible 
algebraically special types, the Petrov type cannot be 
III because of = ^4 = 0, and the Petrov II and N 
cases also have the non-physical equation of state p = 
—p. Finally, due to a theorem by Collinson [l2j], the 
conformally flat case is also excluded. 

Hence, the only algebraic special solutions of physical 
interest one might hope to find are of Petrov type D. 
Up to second order in the rotational parameter equation 
(f2"Tj) gives only one real condition. By substitution of the 
zeroth and first order field equations and ^ into (|2~lj) 
the Petrov type D condition gives the relation pd| 



„dB 
rB — 

dr 



1 — B \ (h 2 - ma) = 



d_ 

dr 



B 2 uj 



(22) 



C. Field equations for the interior region 

Note that mo and ho do not appear in equations ([3]), 
©, ©, ©, JSJ) and JUJ and hence this subsystem for 
A, B, u>, rri2, k% and h 2 decouples. Notice that these 
equations contain m 2 only algebraically. In section [V] 
the system will be reformulated as a coupled system of 
six first order ordinary differential equations. Due to the 
requirement of a regular centre the solutions to this sub- 
system will only depend on four constants of integration. 



D. Vacuum metric 

In the exterior vacuum region we will use a frame 
adapted to the asymptotically non-rotating observer. 
Solving the field equations detailed in Section III Al by 
imposing p = p = the metric functions for the vacuum 
region are given as follows [l], [H| 



.4- 



B 2 = 1 - 2M/r 



2aM 



(23) 



(24) 



ho = 



1 



r-2M 



a 2 M 2 r 



2M 



mo 



2M-r 



1 (o?M 2 



C-2 



C2 



3c x r (2M - r) log I 1 - — ] + a 2 ^- (M + r) 



M 



+2ci— (3r 2 - 6Mr- 2M 



r j r* 

r - M 



2M -r 



2M\ 

,2 OA,f2\ 



r 2 qi 



3d(r 2 - 2 AT) log 1 



M 



2M 



a 2 ^r(2M + r) 



-2 Cl — {2M 2 - 3Mr - 3r 2 ) + [2M 2 - r z ) q x 



m 2 



6a — : h 2 

r 4 



(25) 



In this approximation, the slowly rotating solution is 
characterized by the mass M, the first order small rota- 
tion parameter a, and the second order small constants 
Ci, C2 and q\. When q\ takes the value zero the metric 
is known to be the general asymptotically flat stationary 
and axisymmetric vacuum metric to second order (see 
e.g. [l||). It can be checked by plugging the vacuum 
quantities into the Petrov type D condition |22|) that the 
solution is of Petrov type D only if both ci and q\ are 
zero. The metric is then equivalent to the Kerr metric to 
second order with mass M — > M — c 2 - 

When qi ^ the metric cannot be asymptotically flat. 
It is important to keep in mind, however, that without 
the inclusion of this constant the matching conditions on 
the zero pressure surface are overdetermined in general 



III. MATCHING 

The matching of the fluid ball to a suitable exterior 
vacuum region happens at the zero pressure surface. Be- 
fore matching these two spacetime regions it is informa- 
tive to investigate first the structure of the constant pres- 
sure surfaces. 



A. The constant pressure surfaces 

The pressure in the rotating fluid configuration is given 
by the function p(r, 9) = po{r) +p 2 (r, 6). The surfaces of 
constant pressure, Sf, may be labelled by the function f 
defined by the relation 



p(r, 9) = p (r) = po{r) + 8p{r) , 



(26) 



where po is the corresponding pressure for the non- 
rotating configuration and 8p(f) is a second order small 
shift of the pressure that changes monotonously from the 
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centre where it takes the value £>2o(0) to the zero pressure 
surface where it becomes zero. The value of the central 
pressure to second order follows by assuming regularity 
at the centre and by making use of the field equations 
(fT0|) and (fl9|) . along with the relation for p 2 o (fl3|) (see 
section HV Bp . It turns out that it will depend, among 
others, on one freely specifiable constant corresponding 
to a spherically symmetric perturbation that produces a 
second order small change of the central pressure (and 
density). If we choose to consider only rotational per- 
turbations with P2o(0) = 0, Sp(f) may be chosen to be 
identically zero. 

The radial displacement £ is defined by 

t = r + £ . 

To second order one has 

p{r,8) = p (r) +p 2 (r,9) 

= p (r) +P2o(r) + P22(r)P2{cos9) 



f-\ , t d P° 
= p {f) + 5p(f) , 



Pw(r) +P22(r)P 2 {cos9) 



implying that £ possesses the form £ = £o + £2-P2(cos#), 
where £o and £2 are given as 



£0 = - [P2o(r) - Sp(r)} /(dpo/dr\ f ) 



and 



giving 



h = 2nf 1 + fc 



£0 



ha 
2 



2f 



The length of the curve 7, yielded by the intersection of 
Sf and a plane including the axis of rotational symmetry 
is given as 



£0 
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l 2 = 2irf ( 1 + fc + — + -f + ^ 
r 4 4r 



where we have used the relation 

dl 2 = (1 + 2k)r 2 d9 2 = 
(l + 2fc + 2fc 2 P 2 (cos(9))(f 

along with the fact that the term obtained by substi- 
tuting dr = —3^2 cos 9 sin 9d9 into the line element (JTJ) 
is of fourth order and is hence dropped. The constant 
pressure surfaces are oblate iff l\ > l 2 , i.e. whenever 



-Zo+&P2(cos9)) 2 d9 2 , 



k 2 + ^<0. 



(30) 



Another way of determining the oblateness of the con- 
stant pressure surfaces is possible by comparing the radial 
distance from the origin to the curve 7 and the analogous 
distance to second order in the eccentricity parameter e, 
r = a(l — ie 2 cos 2 9), for an ellipse in M 2 with semi- major 
axis a. The curve 7 is then found to be an ellipse up to 
second order with e 2 given as 



-P22(r)/{dp /dr\f) 



(27) 



Note that there will be a certain arbitrariness in £0, the 
average shift of the radius, unless 8p(f) is specified. At 
the origin £0 = and on the zero pressure surface T\ the 
expressions §5$ and (flU)) together with Sp(ri) = gives 



Co 



12rvl 2 — - - 12to — (rA 2 
dr dr 



1 / 'dv\ 
dr J 



(28) 



D 



m 2 ,_ 

dr 

B 



For infinitesimally small values of r (|3T|) reduces to 



^^ 2 
m 2 + — 

r 



according to which the constant pressure surfaces are 
oblate iff 



m 2 + — < 0. 

r 



(32) 



for Co- If we choose Sp(f) = this expression, with r\ 
substituted with F, holds for any f in the interval [0, r{\. 
From (|5|) and (I14p £, 2 (f) is given by 



& = - 



(3A 2 /i 2 + r 2 ^ 2 ) 



(29) 



The circumference of the intersection of a constant 
pressure surfaces Sf and the equatorial plane 9 = ir/2, 
which is in fact a circle, is obtained from 



Notice that this inequality, along with the relations (|6Tj) 
and (|62|) in section IIVBI does also justify that the two 
different characterization of oblateness are compatible. 

Note, finally, that dr is a form field orthogonal to the 
constant pressure surfaces Sf. Up to second order the 
corresponding normalized field is 



n a = {0,{l + m)/B,3&sm9cos9/B,0) 



B. The matching 



(33) 



dl 2 = (1 + 2fc)r 2 sin 2 



= (l + 2fc -fc 2 )(r + Co-y ) V 



In this section we match the interior rotating fluid so- 
lution to an exterior vacuum region at the zero pressure 
surface Sf= ri . 
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In the vacuum exterior region, suitable hypersurfaces 
for matching are determined by the condition 15] 



si 2 - 



g a = l-C (34) 
To second order such a 



g vv + 2flg vt 4 

where SI and C are constants 
surface can be given as 

r = n 4 x = r \ + Xo + X2^2(cos 6) , (35) 

where xo and Xi are second order small constants. The 
unit normal to this surface is given by 

nW = (0, {l + m^)/B^ v \?,X2S^ecose/B^\Q), (36) 

where the uppercase index W here and after refers to 
vacuum quantities. 

To adjust the coordinates in the two regions we ap- 
ply a rigid rotation in the interior by the transformation 
ip — > ip + Sit . Also, we can rescale the interior time co- 
ordinate by a constant c 4 while matching the spherical 
basis solutions, and then later by a second order small 
constant when doing the matching of the corresponding 
rotating configuration t — > 04(1 + c 3 )t. We do not have 
such a freedom in choosing the time coordinate and ap- 
plying rotation in the exterior region since we want a co- 
ordinate system adapted to asymptotically non-rotating 
stationary observers. 

Together with the values of the other parameters, 
the location of the matching surface Sr 1 is determined 
by the Darmois-Israel conditions [1, [l(|. In particular, 
these conditions pick out the zero pressure surface as the 
matching surface. 

The Darmois-Israel conditions require that the induced 
metrics agree on the matching surface S ri 

ds 2 \s ri = ds 2 {v) \ Sri , (37) 

as well as, the induced second fundamental forms 

K {v) 



K\ s . 

where K is defined as 



\S ri : 



(38) 



K = K ab dx a dx h = h a c h b d n^ c . d )dx a dx 



with 



n a n b + 8 a 



being the projection operator onto the hypersurface or- 
thogonal to the normal vector n a . 
Writing g ab and K ab as 



(o) , 

9ab = g K ab '-\ 

K h - K {0) 



o (1) 

y a b 



o (2) 

y a b 



ah ' ab ' 

where the superscripts W,W ( 2 ) indicate zeroth, first and 
second order terms respectively, one obtains to second 
order 



lab{r) 



9 l >) 



■9$<r)+g%\r) 



- „(°) 



gTin+O+^^+g^in) 



$0i) + 



Or 



(39) 



and similarly for K ab , on the matching surface given by 
r = r% + £• An analogous result holds in the outer region. 
From (|3"7]l . p8| and dSSJ) we then obtain the following 
equations to order by order. 
Zeroth order: 

= c 4 A , SW = b , A {v } r = c 4 A, r (40) 

First order: 

oj {v) = c 4 {uj - SI) , (jW = C A UJ. r (41) 

Second order: 

h[ v) = h 2 , k { 2 v) = k 2 , h { v) =h a +c 3 (42) 

Ci A{h[ v \ - ho, r ) 4 Ho(A {v) . rr - c 4 A, rr ) = (43) 



C 4 rA t r(k 2 r — &2,r) = 



B(rni v) -mo)=to(BW-B, r ) 



(44) 



(45) 



B{m 2 v) - m 2 ) = UB M r - B >r ) + rB{k 2 v) r - k 2<r ) (46) 



XO — £o , X2 — £.2 



(47) 



All quantities here are evaluated at r = r\, i.e. at the 
zeroth order radius. From the zeroth order equations we 
solve for M and C4 as 



M=~n(l-B 2 ) , c 4 = | 



(48) 



The third equation is equivalent to the zero pressure con- 
dition and also gives the radius n implicitly from 



.4 



(49) 



2b2 ±a 

ar 

To first order we solve for a and SI 

Br? dio „ r\ duj , 

From six of the nine second order equations we can solve 
for ci, c 2 , c 3 , qi, xo and X2 as 
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B 



\B 2 {B 4 - 3) (jgy + 36A 2 h 2 {l - B 4 ) + 72A 2 B 2 {h 2 + k 2 ) 



ci 



9{B 2 - lfA 2 r 



(51) 



£o / 9 d-B\ , B 2 { du 



C3 



36A 2 V dr 



C2 



nB 2 (i-B 2 ; 



- /l 



(52) 
(53) 



9rM2(g2 _ 1)6 [I8fc 2 ^(^ 4 - 1)(^ 4 - 8B 2 + 1) 

+21QA 2 B 2 hxB (2B 2 (h 2 + k 2 ) + (1 - B 4 )h 2 ) + 
36h 2 A 2 B 2 (B 2 - 1){B 4 + B 2 - 8) + 

r 4 2 b 2 ((B 2 - 1)(2 + 11B 2 - 75 4 )+ 

6S 2 (B 4 - 3) lnS)] 



Xo = £o and %2 = 6 



(54) 
(55) 



where £o and £2 are obtained from ([28)1 and (|29|) . The 
remaining three equations turn out to be identically sat- 
isfied due to the other matching equations and the field 
equations. Note that we did not assume the Petrov type 
D condition (|22[) when calculating the matching condi- 
tions. Hence an appropriate matching can be done, i.e. 
the vacuum metric in section III Dl is general enough for 
describing the exterior of any axisymmetric rigidly rotat- 
ing perfect fluid ball up to second order. 



IV. NUMERICAL INTEGRATION 

In this section we provide a reformulation of the field 
equations which is more suitable for numerical integra- 
tion. By doing this we can get higher precision at the 
origin where apparent singularities arise, moreover the 
freely specifiable constants are identified more easily this 
way. 



A. Integrating the zeroth order field equation 

In order to simplify ([3]) it is convenient to redefine the 
functions A, h, m and k in terms of the function v as 

A=e u , h = he- 2v , m = me- 2v , k = ke- 2v (56) 

giving 



1 



ds 2 = e 2u {l + 2he- 2u )dt 2 - (l + 2me- 2u )A;dr 2 

B 1 

-(1 + 2ke- 2u )r 2 \d9 2 + sin 2 9 {dip - ojdtf 



The equations simplify considerably due to the fact that 
only the derivative of v will appear. Hence we introduce 
the function z by 



z dv 
B ~ T dr 



1 . 



(57) 



Then the zeroth order equation becomes first order 
in z and algebraic in B [161 ] 



„ dZ rj 

Br— + 2B 2 

dr 



ABz +1 = 



(58) 



furthermore, the pressure of the non-rotating configura- 
tion (El) takes the form 



1 



Po = — (2Bz — B 2 - l) 



(59) 



B. Series expansion around a regular centre 

For sufficiently regular configurations close to the cen- 
tre the metric coefficients can be given as power series 
in r. Assuming that the central pressure and density are 
finite it follows that B(0) — z(0) = 1. The assumption of 
smoothness of the configurations at the symmetry centre, 
in the spacetime sense, implies that the odd coefficients 
in the expansions of the basic variables are zero. Al- 
though the smoothness implies the vanishing of the odd 
coefficients without the use of the field equations, the 
requirement of smoothness of central density and pres- 
sure, together with the field equations, also implies the 
smoothness of the metric functions. 

The vanishing of the odd coefficients can be shown 
in the generic case by plugging power series expansions 
of the dependent variables into the field equations. In 
these cases the smoothness of the density and pressure 
results from these considerations. If odd powers of r are 
included in the expansion of the metric variables then 
it can be shown that the field equations and the Petrov 
D condition, together with the assumption of finite and 
positive central pressure and density, imply the vanishing 
of the coefficients of all odd terms. However, odd terms 
may appear in special cases, and then the metric ceases 
to be smooth at the origin. For example, consider the 
spherically symmetric case when the only field equations 
is given by (j58[) . If r 3 terms are included in the expansion 
of B then one also obtains a non smooth density gradient 
^2.| r= o ^ 0. However, since the field equation implies 

that the pressure gradient ^-| r= o is always zero, the 
squared speed of sound is vanishing at the origin, i.e. 
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v 2 = -j^\ r= o = 0. It is of interest to find the minimal 
requirements needed to guarantee that the solution is 
regular to any order. A conjecture is that this is the 
case when the central pressure and density are finite and 
^£| r=0 ^0 17]. If the Petrov D condition is used instead 
of an equation of state to specify the configuration then 
it can be shown that ^| r= o = can occur only if either 
the central pressure or density is negative. 

Hence, assuming a smooth centre in the spacetimc 
sense, the odd powers will be omitted hereafter. Plugging 
the expressions 



D 


= i + 


hr 2 + b 2 r i + . 




z 


= i + 


z\r 2 + z 2 r A + . 






= U - 


h wir 2 + u) 2 r A - 




h.2 


= 


+ 4V + hi 2 




fh 2 


(0) , (1) 2 

= m 2 + fn 2 r + m 


(2), 
2 


h 




+ k ( 2 1] r 2 + k {2] 


• 4 

r . 



This shows that for realistic configurations bi < and 
z\ > 0, consequently the z\ — b\ term in the denominator 

( 2) 

of h 2 is nonvanishing. Also, the existence of a local 
maximum of the density at the center implies hi > 0. 
The pressure always has a local maximum at r = if the 
central values are positive. 

Assuming hi = implies that the higher coefficients 
in the expansion of h 2 , m 2 and k 2 are zero. We conjec- 
ture that h 2 = m 2 — k 2 = 0, which is also supported 
by a numerical calculation. From this it follows that 
uj = ujq =constant and that A = B = yl + Cr 2 . But 
this simply gives the de Sitter or anti de Sitter solutions, 
depending on the sign of the integration constant C, in 
a rotating frame. 

Plugging the expansions of the original (no tilde) sec- 
ond order spherical perturbation quantities 



- 1,(0) 



h = % 



(60) 



m 



™(!) 2 
Trig r 



,( 2 U 



™( 2 ) 4 

777 r 



(67) 



into the field equations justifies then that all coefficients 
can be given in terms of b\, z\, loq and hi = h 2 . To 
zcroth order one obtains 



4o) = m (o) =A .(o ) = 0i 



(61) 



then to second order 



-hi, wx = -ujq(zx - 36i) , (62) 
5 



and finally to fourth order 

b\ 3u> 2 . , 2 z\ 2 

b 2 = — -H — (Z\— 3oi) , Z9 = Oizi - — of 

2 2 50/7i V 17 ' 2 2 1 



,( 2 ) 



hi(3zi - 136i) 
14 

w 2 (zi - 36i) [2/ii(222;i - 316i) + 3wg(zi - 36i)] 



210/ii(2i - 61) 
fcf = f (* 1 -36 1 ) + ^(61 -21)- />< 2 \ 



/??. 



(2) 



3 

w (zi - 36i) 



,(2) 



^'2 



70/7i 



[/ii(zi — 336i) — 3wg(zi — 36i) 



(63) 



The expansion of the density and pressure of the non- 
rotating configuration can be written as 



Po = Poc ~ 20fi~ ( p ° c + P0cfr 2 + 0(r 4 ) 
1 



Po = Poc - (3po c + 4p 0c p 0c + p„ c ) r 2 + 0(r 4 )(65) 

where the central density and pressure are given by 

p Qc = -661 , p Qc = 2zi . (66) 



together with (|60p into the two remaining equations (jTT 

and (TT5]) gives that two constants, e.g. h^ and mj 1 ' are 
freely specifiable, whereas the other coefficients can be 
expressed in terms of these two. Of these only mj 1 ' is 

essential since the constant /ig can be absorbed by a 
second order rescaling of the time coordinate. To second 
order one gets 



(0) 



(1) _ 5771^/11(21-61) 777 " 



777 ' = and h Q = 



(i) 



2wg(«i-36i) 



(68) 



unless 21 = 36i corresponding to poc = ~Poc- From this 
result it follows that p 2 o, as given by (fTB"]) . at the origin 
is 



10777 n 

p 20 (r = 0) = 



(i) 



hi (61 



J 2 (36i 



zi) 



(69) 



Hence the central pressure will be unchanged if = 0. 
However, the higher order coefficients in the expansions 
(|67p are still nonzero, i.e. the expansions of h and mo 
starts with r 4 terms. We note that, in general, p 2 o(r — 
0) = 6m X) and p 20 (r = 0) = 4h£ ) - 2m 1) . 

Purely spherically symmetric perturbations, corre- 
sponding to a small change of central pressure but un- 
changed equation of state, of the type (|67|) can be ob- 
tained. This kind of perturbations are possible even when 
there is no rotation at all. By choosing luq = hi = and 
/ a spherically symmetric perturbation, with a 



(64) second order shift of the central pressure given by 



6m^j (61 - zi) (3&i - zi) 
5bj + 10b 2 



V2o{r = 0) 



(70) 



is produced. This expression remains valid in general, 
even when the Petrov D condition is not assumed. 
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C. System of differential equations 

Motivated by the results of the previous section it is 
advantageous to define the new dependent variables P, £, 
y, h, k and m through 



B 

z 



l + r 2 /3 
l + r 2 C 
2ry 



h 2 = r 2 h 



(71) 



k 2 = r 2 (r 2 k-h) 
r 2 {r 2 m — h) . 



m 2 



The closed subsystem of equations (J3j) , ((6|), {7}, 
J|) and (|22|) then takes the form 



dr 



r(2/3(/?-C) + C 2 ) 



/3r 2 + 1 



dp _ 2w 2 (C-3/3) 
dr 



3m 



2u 2 r (f3r 2 + 1) 



— \P T 



r_p (c - m 

{Pr 2 + 1) 



du 

— = 2yr 
dr 



dy 
dr 



dh 
dr 



dk 
dr 



~2 

ry 



3m — 5ujy 



ruj (Pr 2 + 1) 



(2u) 



yr 2 ) 



ry(3m- 10/3w 2 (2 + /3r 2 )) 
2iu 2 (f3r 2 + l) 2 
m(3h + uj 2 ) - Au 2 {(h + %- 20k) 



~2 

ry 



2rcu 2 (( - P){pr 2 + 1) 
2uj 2 + 3h\ [Pr 2 + 1) 
3lo 2 (C-/5) 
r 2 m(C - P) 2 - ph(2( - 3/3) 



4fc 



(C-/3)(/3r2 
\-2h(P — Q 



1) 

r (2fc (C - 3/3) + mC! 



r (/3r 2 + 1) 



/3/ 



1 



while m can be solved for algebraically as 



(72) 



2 -u 2 (C-3P)+ 2 -r 2 f(Pr 2 + 1) 2 + 2 -^f- 



y(pr 2 + l) 2 [2w(C - 3/3) - y (l + /3r 4 (C - /3) + 2/3r 2 )] 



2r 2 w 2 



/3(C - 2/3) ( 



3h 



r 2 (Pr 2 + 1) [w 2 (C - 2/3) - y(/3r 2 + \){r 2 y + 2u))] +3h - Pr 2 u 2 
r 2 uj 2 (C-3P)) -c 2 (C-3/3) 2 



r 2 (/3r 2 + 1) [w 2 (C - 2/3) - j/(/3r 2 + l)(r 2 y + 2w)] + 3h - Pr 2 w 2 



Boundary conditions at r = are given as 

/3(0) = 6i, C(0)=2i, ^(0)=c^o, 
2 

y(0) = wi = -u} (zi - 3bi),h(0) = fti, 
5 

fc(0) = k?+h^ = ^iz 1 -3b 1 )+ } ^ib l ~ Zl ) . 
The relation 

m(0) = mf 5 + 4 2) = ^o(^i ~ 3& i) 

is then satisfied identically. As we have seen there are 
four freely specifiable constants: b%, z\, u>o and hi- 

The system of equations (f?2"j) possesses two types of 
scale invariances. The first one is associated with the 
rescaling of the r-coordinate, r — > ar, under which trans- 
formation the dependent variables scale as 

/3,C,m, k,y^ 2,-0,-0, u,h->Lj,h. (73) 

or a z a z a z a z 



There is also a rescaling associated to the rescaling of the 
rotational parameter loq, following the rule luq — > 70^0 , 
which induces the transformation 

^,y-^l^,iy, h,rh,k -> -f 2 h,j 2 m,j 2 k, 
P,(,r^PX,r ■ (74) 

It is interesting that a combination of the above two 
rescalings with 7 = 1/0? yields a similarity transforma- 
tion of the investigated system. 

Due to these scale invariances of the equations two of 
the constants, e.g., ojq and b\, can be fixed. All other 
configurations can be obtained by rescaling. Note also 
that b\ and Z\ can be expressed in terms of the zeroth 
order central density and pressure as 

h 1 1 

»1 = - ~P0c, Zl = 77POc- 
D 2 

Some of the above equations contain terms of the 
type " F/r", thereby they are apparently singular at 
the origin, so we collected them to the beginning of 
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the right hand sides. It can be checked case by 
case that all of the corresponding numerators vanish 
at the origin. Nevertheless, in determining the val- 
ues of the corresponding ratios numerically it turned 
out to be advantageous to use as the fundamental vari- 
ables the differences (a, /3a, ^a, Va, h&, wa, &a between 
the variables (3,oj,y,h,m,k and their exact values 
Co, A) j wo, ?/0j ^0) ko at the origin. Due to the can- 
cellation of the terms involving the exact values at the 
origin what remains from the numerators will be propor- 
tional to r 2 . 



NUMERICAL SOLUTIONS 



The system (|72l) . when rewritten in terms of the 

variables £a, /3a, Wi, VA> h&, wa, &A, was solved using 
fourth order Runge-Kutta. A check of the convergence 
factor 



fn f2n 
f2n — fin 



where n is the number of points in a given r-interval, 
was performed for various quantities /. We found that 
the errors decreased according to the expectations, i.e. 
the value of C n was found to be close to 16. 

We also checked whether the scale invariance proper- 
ties of the field equations were reproduced properly by 
our numerical code. When we rescaled the freely speci- 
fiable boundary data according to (|73[) or ([7^]) then all 
the dependent variables scaled in the appropriate way. 

In scanning the four dimensional parameter space, due 
to the scaling invariances, without loss of generality, we 
fixed bi = — 1 and loq = 0.1 while we varied the central 
values z\ and h\. Notice that to have positive central 
densities and pressures the relations b\ < and z\ > 
also have to be satisfied. The integrations were carried 
out until the zero pressure surface was reached. 



A. Check of Wahlquist 

The code was checked for the Wahlquist solution [l8| . 
which is of Petrov type D. To second order it is given by 



ds 2 = f (l + 2h)dt 2 - 2^-^dx 2 - 
Mo^/o 

sin 2 x(l + 2k) \d9 2 + sin 2 6 {dip - ujdtf (75) 



Mo^ 



with 



and 



/o = H — o 0- _ ^ cot a;) 



2 sin 2 x 



(1 — x cot x) . 



(76) 



(77) 



The transformation to the Hartle variables used in ([T]) is 
given by 



• since 



For the functions h, m and k see [13| . For the Wahlquist 
solution the four starting values are given by 



Mo 
12 



and 



(1-3k 2 ), zi 



hi 



fJ-o 



(1-K 2 ), UJ 



60 



(78) 



in terms of the three integration constants ^o, K and tq. 
Solving for hi gives 



hi 



(zi - 3bi)u;l 
5(6i - zi) 



(79) 



Note that hi < for positive central density and pres- 
sure, implying together with equation (|64[) the well 
known fact that the density of the Wahlquist solution 
has a minimum at the centre. 

In figure [T] the relative error between the analytical 
solution for the the rotational function u>, as given by 
(|77| . and the numerical solution is plotted for various 
resolutions. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



FIG. 1: The relative error |aj numcri cai — ^analytic | /^analytic is 
shown for the resolutions Ar = 0.01, Ar = 0.005 and Ar = 
0.0025. Starting values are zi = — bi = 1 and ujq — 0.1 giving 
hi = -0.004. 



As it was already found [l3|, [l9[ the shape of the 
Wahlquist fluid ball is always prolate which is also in ac- 
cordance with the positive sign of the quantity ^-(-^Ai- 
The quantity qi is not zero for the Wahlquist solution 
which is equivalent to that it cannot be matched to 
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an asymptotically flat exterior solution to second order 
[l3l [l9l . [20| . The value of c\ was also found to be nega- 
tive for all the tested Wahlquist configurations. This is 
also verified by analysing the analytical expression for c\ 
given in [19( . It would be interesting to know whether the 
sign of this quantity is related to the shape in general. 



q, =0 



B. Asymptotically flat solutions 



A solution to the field equations (|72|) is asymptotically 
flat iff qx = 0. In Figure the points of the dashed 
curve represent configurations with value q\ = in the 
Zi/ii-plane, while in Figure [3] a section of the same curve 
for small Zi, corresponding to small central pressures, 
is given. Naturally, the asymptotically flat solutions, 
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0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 
hi 

FIG. 2: Along the solid curve the density and pressure become 
zero for the same value of the radius r = ri . Below this curve 
those configurations can be found which can be matched to an 
exterior vacuum region. The dashed curve represents those 
configurations in the zi/ii-plane for which the exterior vacuum 
region is asymptotically flat. 

represented by points belonging to these curves, always 
have finite radii, but further increasing z±, correspond- 
ing to the increase of the central pressure, we get into a 
region where the density becomes negative before a zero 
pressure surface is reached. The limiting curve where the 
pressure and density become zero at the same radius is 
shown in Figures [2] and [3] by the solid lines. Due to the 
negative density the pressure ceases to be a monotonic 
function of r for configurations above the p ri = curve. 

Since for the asymptotically flat solutions hi may be 
seen as a function of zi, we see that h 2 , 1112 and k 2 for 
these solutions are determined by 61, z\ and ljq. The 
analysis done in section lTVBI for the remaining field equa- 
tions showed that essentially one constant of integration 
is freely specifiable for the functions h a and m . This 
implies then that the second order configuration is com- 
pletely determined by the zeroth order spherical config- 



0.5 - 

— ' 1 — ^ 1 1 1 1 1 1 

0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 

FIG. 3: The same curves as in Fig. [2] are shown for small z\, 
corresponding to low central pressures. 



uration, a second order spherically symmetric perturba- 
tion and the magnitude of the rotation. 



C. Are there Kerr-like solutions? 

To our knowledge the only known source for the Kerr 
metric is the thin rotating disk of dust with a = m found 
by Neugebauer and Meinel [2l[ . It is tempting to investi- 
gate whether there can exist a fluid ball belonging to the 
class investigated in this paper that could be matched 
to the Kerr solution to second order. To settle this is- 
sue note first that the metric of the exterior region be- 
comes the Kerr metric with mass parameter M — c 2 iff 
qi = c\ = 0. However, the numerical runs indicated that 
either c\ > for hi > or c\ < for hi < occurs, 
in general, i.e. the desired matching seems not to be 
supported. Note that these numerical findings are in ac- 
cordance with some earlier results, see e.g. UH2|j telling 
that typically the exterior metric deviates from the Kerr 
metric due to the ellipsoidal shape of the rotating fluid 
ball. To this end it is illuminating to consider an ex- 
pansion of the exterior metric for large r which gives the 
following leading terms of <?oo (with qi = 0) 

1 2M(l-ff) 2MP 2 {cos9)(a 2 + f M 4 Cl ) 
.9oo — 1 1 5 , 

(80) 

i.e., the associated quadrupole moment reads as (cf., e.g. 

M) 

Qn = Q22 = -Q33/2 = -2M (V + yM 4 Cl j 

in an asymptotically Cartesian system with the 3-axis 
along the axis of rotation. In Figure 2] the value of ci 
as function of the central pressure p 0c — 2z% along the 
qi = curve is shown. 
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10 

Pa, 



FIG. 4: The constant ci as a function of zeroth order central 
pressure, poc, along the curve qi = 0. 



D. Some asymptotically flat solutions with 
reasonable equation of state 

In this subsection we present some properties, like 
equation of state and speed of sound, for some of 
the physically interesting inner fluid ball configurations 
which can be matched to a suitable asymptotically flat 
exterior vacuum region up to second order, i.e. those 
solutions for which qi vanishes. According to TableU 
the value of hi has to be smaller than around 0.012 
for these solutions to have subluminal speed of sound, 
Vg = dp/ dp < 1. For all configurations in this interval 
the speed of sound also increases when approaching the 
centre, i.e. the fluid becomes stiffer as would be expected 
on physical grounds. 

In Figures [5] and [6] the zeroth order pressure and den- 
sity, po and po , are shown as functions of r for some con- 
figurations with central pressure po c between 0.218 and 
4.366, while in Figure [7] and [8] the equation of state, i.e. 
p as function of p, and the square of the speed of sound, 
v s = 3^) as a function of r are depicted, respectively, for 
the same family of solutions. 

Notice that the value of v 2 s at the surface of the fluid 
ball seems to be universal, i.e. apparently independent of 
the values of the free parameters at the centre. Unfortu- 
nately, we could not find an analytic argument support- 
ing this observation, nevertheless, it is interesting on its 
own right that such a lower limit may exist for the class 
of the investigated fluid balls. 

In general, the equations of state cannot be polytropic 
since the density does not tend to zero while approach- 
ing the matching surface. However, the equations of state 
can be approximated with a polytropic one close to the 
centre. The two last columns in TableQ] provides the cen- 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



FIG. 5: The zeroth order pressure po is shown as function of 
r for various central pressures along the curve qi — 0. From 
top to the bottom the central pressures are given by 4.3666, 
2.6978, 1.7762, 1.4872, 1.21, 0.944, 0.6896, 0.448 and 0.2184 
respectively, whereas the central density is poc = 6 = — 6&i. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



FIG. 6: The zeroth order density po is shown as function of 
r for various central pressures along the curve qi = 0. From 
top to the bottom the central pressures are as given in figure 

m 

tral values of the adiabatic index 

ndp p + pdp 

k=-— = — , (81) 

p dn p dp 

where n is the baryon number density, and the Newtonian 
adiabatic index 

- P d P (QO\ 
K N = -—, 82) 

pdp 
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Using the Taylor expansion around the centre, as it was 
made in subsection llV Bl one obtains in the limit f — > 



FIG. 7: The equation of state p — p(p) is shown for the same 
configurations as before. 



S. 1.2- 




FIG. 8: The square of the speed of sound, = -f-, is shown 
for the same configurations as on the previous figures. 



that approximates K for low pressures, respectively. Un- 
fortunately, the value of k is not in the preferred range 
4/3 — 5/3 that is considered to be physically acceptable 
in the case of compact neutron stars or white dwarfs. 

As we have seen in subsection llll A[ the surfaces of 
constant pressure, S r , are determined by the relation r = 
f + £q + £ 2 -P2(cos0), where £ 2 is given by ([^5)1 . These 
surfaces are oblate iff fca(f) + £2/^ < 0. In terms of the 
functions v, (3, C? h, k and u> the expression k2{f) + £2/^ 
can be written as 



j / \ j. 1 2v it 2? {3h + lu 2 )(1 + f3r 2 ) 

k 2 {r)+£, 2 /r = e~ 2v I r 4 fc -r 2 h- v 



lim k 2 (r)+£ 2 /r 
1 — >o 



a "<°>(wo+3/u)/(*i-&i) . (84) 



Note that e u ^ may be fixed to 1 since its value only 
corresponds to a rescaling of the time-coordinate. Since 
Z\ > and b\ < for realistic configurations we see 
that close to the centre the surfaces of constant pres- 
sure are oblate iff hi > —uJq/3. In TableQ]the values of 
a = (k 2 (f) + &/f) |f=n and b = (fc 2 (f) +&/?) |f=o are 
also indicated for several configurations. In virtue of the 
negative signs of these parameters the constant pressure 
surfaces are all oblate for these configurations. 

In TableU the central pressure pa c — 1z\ , the radius 
of the zero pressure surface r\, the shape at the zero 
pressure surface, the shape close to centre, the value of Ci, 
the maximal speed of sound v 2 sl the zeroth order density 
at the zero pressure surface and the adiabatic index are 
given for a sequence of solutions with q\ = 0. 



3(C-j8) 



(83) 



VI. CONCLUSIONS 



The most important finding of this paper is that a sub- 
class of slowly rotating perfect fluid balls of Petrov type D 
can be matched to asymptotically flat vacuum spacetimes 
and also that in general slowly rotating perfect fluid balls 
can be matched to non-asymptotically flat vacuum exte- 
riors determined by equations (j23)) . (|24|) and (|25| . Our 
numerical results support the conclusion that neither of 
the Petrov type D inner fluid solutions can be matched to 
second order to the Kerr metric, which is in accordance 
with the general expectation that the ellipsoidal shape 
of the rotating fluid ball produces an extra contribution 
to the quadrupole moment which should also be present 
in the corresponding quadrupole moment of the external 
field 0,(13. It was also found that there is a range in pa- 
rameter space for which the value of the central pressure 
is relatively low and the speed of sound is also sublu- 
minal. The equation of state was also determined for 
various solutions belonging to the investigated class. It 
is clear that the equation of state cannot be polytropic 
since, in general, the energy density does not vanish at 
the zero pressure surface. Nevertheless, the equation of 
state can be approximated close to the centre by a poly- 
tropic one. Unfortunately, the corresponding adiabatic 
index k was found to take values out of the physically 
preferred range. 
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/ll 


poc — 2zi 


ri 


a 


b 


Cl 


2 
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po(ri) 
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KN 


0.006 


0.2184 


0.264 


-0.0086 


-0.00841 


1.13 ■ 10 3 


0.36 


5.3021 


10.3 


9.9 


0.007 


0.4480 


0.367 


-0.0089 


-0.00844 


46.675 


0.44 


4.7244 


6.3 


5.9 


0.008 


0.6896 


0.441 


-0.0091 


-0.00843 


8.165 


0.54 


4.2416 


5.2 


4.7 


0.009 


0.9440 


0.500 


-0.0093 


-0.00839 


2.556 


0.64 


3.8318 


4.7 


4.1 


0.01 


1.210 


0.548 


-0.0095 


-0.00831 


1.105 


0.74 


3.485 


4.4 


3.7 


0.011 


1.4872 


0.587 


-0.0098 


-0.00822 


0.584 


0.85 


3.1898 


4.3 


3.4 


0.012 


1.7762 


0.621 


-0.0100 


-0.00812 


0.352 


0.97 


2.9324 


4.2 


3.3 


0.015 


2.6978 


0.700 


-0.0109 


-0.007805 


0.1259 


1.35 


2.3441 


4.4 


3.0 


0.02 


4.3666 


0.786 


-0.0126 


-0.00733 


0.0514 


2.05 


1.7335 


4.9 


2.8 


0.05 


15.2 


1.029 


-0.0254 


-0.00620 


0.0193 


6.76 


0.66428 


9.4 


2.7 


0.1 


32.8 


1.190 


-0.0449 


-0.00594 


0.0210 


14.95 


0.42164 


17.7 


2.7 


1 


312.8 


1.636 


-0.2297 


-0.00637 


0.0503 


164.58 


0.21398 


168 


3.2 



TABLE I: The central pressure po c = 2zi, the radius of the zero pressure surface n, the shape of zero pressure surface 
(a = + £.2/ri is negative for an oblate configuration), the shape of the constant pressure surfaces close to centre given by 
b = (fe+$2/f)|f=o, the value of the constant ci, the maximal speed of sound Vg, the zeroth order energy density at the matching 
surface and, finally, the central values of the adiabatic indices k = and kn = are shown for some configurations 

which can be matched to an asymptotically flat exterior (with q\ = 0) . Although for all listed configurations the central density 
Poc = 6 = — 66i any central density can be obtained using the rescaling freedom (|73(l . 



